Elastic constant dishomogeneity and Q 2 dependence of the broadening 
of the dynamical structure factor in disordered systems 
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We propose an explanation for the quadratic dependence on the momentum Q, of the broadening 
of the acoustic excitation peak recently found in the study of the dynamic structure factor of many 
real and simulated glasses. We ascribe the observed Q 2 law to the spatial fluctuations of the local 
wavelength of the collective vibrational modes, in turn produced by the dishomegeneity of the inter- 
particle elastic constants. This explanation is analitically shown to hold for 1-dimensional disordered 
chains and satisfatorily numerically tested in both 1 and 3 dimensions. 



PACS numbers : 63.50,+x, 61.43.Fs 

The study of the nature of collective atomic excitations 
in disordered solids at wavelengths A approaching the in- 
terparticle separation a, a/ A ~ 0.1 — 0.5, has received 
renewed interest in the last few years thanks to new ex- 
perimental tools and to improved numerical techniques. 
On the experimental side, the dynamics of the collective 
excitations is often investigated via the dynamic struc- 
ture factor S(Q,uj), i. e. the space-time Fourier trans- 
form of the particle density correlation function. The 
study of S(Q,u>) in the region of mesoscopic exchanged 
momentum (Q = 1 — 10 nnr 1 ) has become recently pos- 
sible in disordered systems thanks to the development of 
the Inelastic X-ray Scattering (IXS) technique jL|, and 
many glasses (2j and liquids || have been studied with 
this technique. Although specific quantitative differences 
exist among different systems, all the investigated glasses 
show some common features that can be summarized as 
follows: i) there exist propagating acoustic-like excita- 
tions up to Qa w 3; ii) the slope of the Q — u> dispersion 
relation in the Q — > limit extrapolates to the macro- 
scopic sound velocity; Hi) the broadening of the excita- 
tion peaks in S(Q, u>) follows a Q 2 law. These general fea- 
tures of S(Q, uj) have been confirmed by numerical calcu- 
lations in different glasses, using both standard Molecular 
Dynamics (MD) simulations 0], and the Normal Mode 
Analysis (NMA) in the harmonic approximation H . 

The Q 2 dependence of the excitations broadening has 
not yet received a theoretical explanation. Such depen- 
dence is the same as predicted by hydrodynamics, but 
this coincidence is only accidental. Indeed, in the exper- 
iments the broadening is found to be temperature inde- 
pendent and the Q 2 law is numerically found also in har- 
monic glass models; these results indicate that the origin 
of this behavior should be found in structural rather than 
dynamical properties, i. e. it should be associated to the 
atomic disorder in the glass and not to dissipative phe- 
nomena like anharmonicity or relaxation processes. For 
this reasons, we concentrate on harmonic systems, leav- 
ing apart all the difficulties due to dynamical processes 



like, for example, anharmonicity. 

In this Letter we suggest that the observed broadening 
is due to the spatial fluctuations of the elastic constants. 
Our explanation is based on analytical and numerical 
results for the simple case of a 1-dimensional (1-d) har- 
monic disordered system (linear chain of like mass joined 
by random springs), and then it is generalized to 3 di- 
mensions. 

The dynamics of 1-d disordered lattices has been thor- 
oughly investigated in the past, and a recent overview on 
this subject can be found in Q. In principle, 1-d systems 
often show very peculiar characteristics; more specifically 
it is well known that, at variance with d>l cases, in 
a (infinite) disordered linear chain all vibrational modes 
are localized and localization by itself contributes to the 
broadening of the inelastic peaks of S(Q,u>). However, 
as we will show, this effect can be disentangled. Our 
model consists of an linear chain of N particles (i — 1..N) 
of mass M placed a distance a apart from each other 
(xi — ia), and joined by next-neighbour springs (Ki) 
randomly chosen from a flat distribution (V(K)) with 
extrema K ± AK (mean fi^ = K and standard deviation 
(Tk = Aif/vo). The dynamics of the system can be ex- 
pressed in terms of its eigenvalues (u) p ) and eigenvectors 
(e p (i) ), p = 1...N being the mode label. The classical 
dynamic structure factor is expressed as: 



S(Q,u) = ^%E(Q,u) 
Ad uj^ 



(1) 



with 



£(Q,^H5>~ P (Q)| 2 (5( W -^), (2) 
p 

e P (Q) = N- 1 ' 2 exp(iQxi)e p {i). 

i 

The eigenvectors and eigenvalues are the solutions of 
the secolar problem, J2j Dije p {j) = u> 2 e p (i) where Z>y is 
the dynamical matrix. The basic properties of S(Q,u>) 
are cointained in E(Q,uS), the squared space Fourier 
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transform of the eigenvectors of modes with frequency 
ui p "close" to u>. A typical eigenvector obtained from the 
diagonalization of the dynamical matrix is reported in 
Fig. 1. As is well known (see for example there ex- 
ists a rather well defined wavelength but, as compared 
with the sinus function expected in the case of an or- 
dered system, three main differences can be noted: i) the 
peak height is not constant, i. e. there exists an envelope 
which is localized in space (see inset a); ii) the wavelength 
A (i. e. the distance between two next nearest nodes) is 
not constant, but is rather a space fluctuating quantity; 
its statistical distribution is reported in the inset b) of 
Fig. 1; finally, Hi) by analysing in detail the eigenvector 
between two successive nodes (see insets c) and d) , some 
deviations from the simple sinus law can be evidenced. In 
principle, all these three characteristics could contribute 
to broade the inelastic peaks in S(Q, u>). In the following 
we will show that the broadening is mainly due to fluc- 
tuations of the local wavelength and (in 1-d systems) to 
localization, while Hi) has minor consequences, and its 
signature can only be found in the low-w tail of S(Q,uj). 

Let us first focus on ii), i. e. on the effect of spatial 
wavelength fluctuations. To show how these fluctuations 
by themselves can produce the Q 2 -broadening, we build 
up a model eigenmode which possess only characteristic 
ii), that is a sequence of alternating positive and negative 
semiperiods of a sinus function of slightly different wave- 
length. If h m (m = I, 2...) are the positions of the nodes, 
we define the local wavelengths A m = h m — h m -2 and 
compute the Fourier transform of the model eigenvec- 
tors, i.e. e(Q) oc Jh"* 2 ^ x e:K P(^Q x ) sin(27r/ X rn (x — 
h m -2))- This is then averaged over different realizations 
of the disordered chain, assuming a gaussian distribution 
of wavelengths centered at A and with variance o\ . The 
calculation of E(Q) = N^ 1 (\e(Q)\ 2 ) is straightforward 
althought long; it produces an E(Q) that near the peak, 
located at Q = 2tt/X, has a quasi- lorentzian shape with 
half width at half maximum (HWHM) T p given by: 



-,(0) 



A 2 



(3) 



The successive step is to find a relationship between a\ 
and the characteristics of disorder. In each segment of 
length A m , containing n m springs, we define a "local" 

sound velocity, v 1 ^ = ay K%{ J ' /M, where Kffi is the ef- 
fective elastic constant obtained by averaging the individ- 
ual spring constants Ki that are found between h m -i and 
h m : (K^Jf)- 1 = Y,ie m Pi K 7 X /Y,izmPi- Here the P t 's 
are weights that take into account the fact that springs 
near the nodes are more effective (highly stretched) than 
those near the antinodes (not stretched) in determin- 
ing K$[* . The validity of the concept of "local veloc- 
ity" can be tested numerically. In fact, coming back to 
the eigenvectors obtained by the solution of the secular 



problem for the disordered linear chain, the actual wave- 
length of the m-th half-period, A m , can be compared 
with that obtained from the local sound velocity, i. e. 
from the relation A^° = 2irv l £f /u>. As an example, in 
Fig. 2 we show the correlation between A and X loc ob- 
tained from the analysis of the eigenvectors of the modes 
at u/y/K/M fa 0.031 of 1,000 disordered chains with 
AK I A=0.3, for different choices of the weights Pj. As 
can be seen, the correlations are satisfactory, indicating 
the validity of the assumptions that i) a local sound ve- 
locity exists, and ii) it is determined by a local spring 
constant averaged over the wavelength. Moreover, we 
observe that the highest correlation is found for weights 
proportional to the square of the local strain. From the 
previous relations it is easily deduced that: 



T ~~ ~ ~ 2 K e fl 



f <Jk--ss _ 1 facr k 

n T77TT — 7:\ — - 



2 V n 



(4) 



where n = (n m ) = A/a and a — (P 2 )/(P) 2 . In the 
following we will choose P proportional to the square of 
the strain and, therefore, a — 3/2. By substituting Eq. || 
into |^, we obtain: 



r (Q) =a^-0 2 



(5) 



which reproduces the experimental Q 2 behavior. Sum- 
ming up, we have shown here that the spatial elastic 
constant fluctuations by themselves produce a broadened 
E(Q,oj) which is quasi-lorentzian in shape with HWHM 

rf oc q 2 . 

In order to check the validity of Eq. || we performed 
a numerical calculation of E(Q,u>). At selected val- 
ues of oj, and for different values of AA", we calculated 
the eigenvectors of 50 different realizations of a disor- 
dered chain composed of 20,000 atoms using the Dyson- 
Schmidt (DS) method ||. Using Eq. || the functions 
E(Q,oj) are then calculted as a function of Q; represen- 
tative examples are reported in Fig. 3 for the indicated 
values of u> and AK/K = 0.6. The E(Q,u>) have been 
fitted to a lorentzian lineshape. The derived HWHM 
r are reported in Fig. 4 together with our prediction 
(Eq. ||), which reproduces correctly the Q 2 behavior but 
is about a factor 2 too small. This discrepancy is due to 
our neglecting of localization effects which, as mentioned 
before, are certainly effective in 1-d systems. In these 
systems the eigenvectors have an exponential envelope in 
the tails, whith a decay length L, and this produces a 
contribution to the linewidth of E(Q,lu) vs Q given by 
r^ Q = 1/L Q. Since localization also gives rise to a 
lorentzian lineshape, we expect a total broadening given 
by r <Q) = r <Q) +r <Q> . The values of L have been numer- 

J Tot L F 

ically computed by fitting the tails of the eigenvectors 
produced by the DS technique; as expected || we found 
L oc lj~ 2 and, therefore, also in this case T Q> oc Q 2 . The 
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full line in Fig. 4 represents the total HWHM F J 3 ' : it 
is in quantitative agreement with the widths of E(Q,u>) 
determined numerically. 

We stress that the effect of localization is peculiar to 
1-d systems and is not relevant to 3-d disordered systems 
H . Therefore, we expect that in 3-d the wavelength fluc- 
tuations induced by the spatial flutuations in the local 
clastic constant account for most of the broadening ob- 
served in experiments. Since usually experiments probe 
S(Q,u) as a function of u> at fixed Q, we rewrite Eq. |5| 
as 



r M -n (Q) a k a rfl 
T F =vT F =av-j-Q 



(6) 



where v is the sound velocity. In order to verify whether 
this equation, derived for 1-d, is still valid for 3-d systems 
we determined numerically the broadening of S(Q,u>) 
for a simple model glass consisting of 32,000 Ar atoms 
interacting via the 12-6 Lennard- Jones (LJ) potential 
{a = 3.405 A, e/K B = 125.2 K); for further details on 
this simulation see |hJ. In the inset of Fig. 5 we re- 
port the r values measured from a lorentzian fit to the 
numerically calculated S(Q,uj) (which are also shown in 
Fig. 5), compared with the prediction of Eq. ^ (full line). 
The values of 07., and a, needed to calculate r ' 1 from 
Eq. ^3], were determined by averaging the interatomic dis- 
tances and the second derivatives of the LJ interatomic 
potential |li"| ] . The good quantitative agreement between 
the numerical and the predicted values of T ( 1 gives us 
confidence on the applicability of Eq. || to real glasses, 
and therefore on the prevalence of the elastic-constant 
fluctuation mechanism in determining the width of in- 
elastic peaks in the S(Q,cv). 

We emphasize that the mechanism proposed for the 
origin of the broadening of the Brillouin peaks in S(Q, u>) 
is completely different from the situation where a mode 
with well defined wavelength exists and its scattering is 
produced by isolated point defects (Rayleigh scattering) , 
which would result in a Q law. 

In conclusion, we propose a plausible explanation for 
the ubiquitous Q 2 law found for the broadening of the 
Brillouin peaks in disordered systems. According to the 
present model this broadening is a consequence of the 
spatial fluctuation of the "effective" (i. e. averaged over 
one wavelength) elastic constant. This fluctuation be- 
comes smaller and smaller as the wavelength (i. e., as the 
number of involved springs) increases and in particular, 
from basic statistics it is oc A -1 / 2 . This result, together 
with the existence of a "local velocity", leads to a Q 2 
dependence of the Brillouin linewidth in harmonic disor- 
dered systems at all Q's. 
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The quantities a, fi^, and <7fc, entering Eq. |(| were deter- 
mined by the equations: 



(0> = 



dr 4-7T r 2 g(r) 0(r) 



dr 4-7T r 2 g(r) 



where 0(r) = r for a, 0(r) — V LJ (r) for fik and 
0(r) = (V^'j(r) - Hk) 2 for a 2 . Here V L J {r) is the LJ po- 
tential function, g(r) is the pair distribution function of 
the system and r a a cut-off radius that takes into account 
that in the theoretical derivation of Eq. ^| we only con- 
sidered nearest neighbours interaction. The results only 
slightly depend on the choice of r , and we fixed it at half 
way between the first and second peak in the g(r). 

FIGURE CAPTIONS 



FIG. 1 Portion of the eigenvector of the mode at o;/-^/ K/M = 
0.571 of a chain with N=1,000 particles and AK/K = 0.6 
reported as a function of the particle coordinate. In the 
insets the three most relevant features of the eigenvectors 
are emphasized: i) Localization (the envelope of the whole 
eigenvector is reported in inset a); ii) Wavelength fluctua- 
tions (the distribution of A m calculated from 100 realization 
of disorder is reported in inset b together with the gaussian 
fit, A/a=10.18, cr\/a=0.827); and Hi) Deviation from sinus 
curve (a blow up of the eigenvector, dots, is reported in the 
inset c) together with the best local sinus approximation, 



3 



full line. In the inset d) it is reported the residual of the 
eigenvector after the subtraction of the full line in the inset 
c).) 

FIG. 2 Countour plot of the numerical joint distribution func- 
tion V(\,\ loc ) obtained by the analysis of the modes at 
w/ \J K/M ss 0.031. The countour plot has been obtained 
averagin 1,000 realization of disordered chains (AK/K = 
0.3) with N =1,000. The analysis has been performed with 
different choices for the weigth: a) P=l; b) P = |Ve p (i)|; 
and c) P = |Ve p (i)| 2 . The P-scale is logarithmic and the 
i-th line indicates the level at 10 — t / 4 Vmax- 

FIG. 3 Examples of E(Q,lj) vs Q at the indicated ui values cal- 
culated with the Dyson-Schmidt method |9] averaged over 
50 realizations of disordered linear chains (AK/K = 0.6) of 
length N = 20, 000. Only the peak region is reported, and 
the full lines are the best lorentzian fits in this region. 

FIG. 4 w-dependence of the hwhm of the lorentzian fit to the 
E(Q,u>) reported in Fig. 3 (full dots) compared with the 
predictions from elastic constant fluctuations including (full 
line) or not including (dashed line) the localization effects. 

FIG. 5 Examples of S(Q,u>) of simulated glassy Argon at the 
indicated Q values. In the inset the hwhm of the peaks of 
the S(Q, lu), calculated by a lorentzian fit in the peak region, 
(full dots) are reported vs Q together with the theoretical 
prediction from Eq. B (full line). 
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